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Abstract 

In this paper we study the quantization of the nonlinear oscillator 
introduced by Mathews and Lakshmanan. This system with position- 
dependent mass allows a natural quantization procedure and is shown to 
display shape invariance. Its energy spectrum is found by factorization. 
The linear harmonic oscillator appears as the A ^ limit of this nonlinear 
oscillator, whose energy spectrum and eigenfunctions are compared to the 
linear ones. 
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1 Introduction 

In a recent paper [1] we have analyzed a classical nonlinear oscillator which is a n- 
dimcnsional generalization of the one-dimensional system introduced previously 
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by Mathews and Lakshmanan [2], [3], as a onc-dimcnsional analogue of some 
models of quantum field theory [4] , [5] . Such mechanical system was described 
by a Lagrangian 



which represents a nonlinear oscillator-like with amplitude dependent frequency 
periodic solutions. Note that this system can also be considered as an oscillator 
with a position-dependent effective mass m = (1 + Xx^)~^ (see e.g. [6] and [7] 
and references therein). 

The two-dimensional generalization studied in [1] was given by the La- 
grangian 



and it was shown to be not only integrable but super-integrable. This suggests 
that the corresponding quantum model should be exactly solvable, although one 
may expect to have some ordering ambiguities because of the a;-dependence of 
the kinetic term. 

We aim to study in this paper the one-dimensional quantum model using the 
well-known techniques of factorization and related operators (see e.g. [8] and [9] 
and references therein). This algebraic technique was started by Schrodinger 
in [10] and [11] and its interest has been increasing since the beginning of Su- 
persymmetric Quantum Mechanics (see [12] for a review). 

In more detail, the plan of the article is as follows: Section 2 is devoted 
to study the simplest a = case in the classical approach [1], both in the La- 
grangian and the Hamiltonian formalism, and an infinitesimal symmetry as well 
as the invariant measure in M under such vector field are determined. Then, 
we proceed to introduce the quantum Hamiltonian describing this position- 
dependent free system. The factorization method for this kind of position- 
dependent mass [13] is developed in Section 3 and the specific example of the 
quantum nonlinear oscillator is studied in Section 4, where we prove that the 
problem has shape invariance. In Section 5 the spectrum of the qiiantum non- 
linear oscillator is found by using the method proposed for such systems by 
Gendenshtei'n in [14, 15]. The last section includes some final comments on the 
relation of this problem with that of the harmonic oscillator in the limit A ^ 0. 



Let us first recall the case of the onc-dimensional "free-particle" motion (in the 
sense that a = 0) characterized by the Lagrangian 




(1) 




2 A-dependent "Free Particle 




(3) 
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As it was remarked in [1], the Lagrangian function Ti(A) is invariant under the 
action of the vector field = (A) given by 

X,(A) = /l + A^^, 

in the sense that we have 

X*(A)(ri(A)) =0, 

where Xl.{X) denotes the natural lift to the phase space MxR (tangent bundle 
in differential geometric terms) of the vector field Xx{X), 

This vector field can be seen as a Killing vector field for the metric g = 
{1 + Xx'^)~^ dx®dx, and generates a one-parameter group of isometries in this 
Riemann space. The natural measure in the real line is not invariant under 
such vector field; instead, the only invariant measures are the multiples of dfj, = 
(l + Ax2)-i/2da;. 

It is important to remark that in the space -L^(M, d/x) of square integrable 
functions in the real line, the adjoint of the differential operator 1 + Xx'^ d/dx 
is precisely the opposite of such operator. 

On the other side, with the momentum defined as usual, 

dL v^. 

P 



dvx 1 + A a;^ ' 

the Legendre transformation {x,Vx) i— > {x,p) leads to a Hamiltonian function 
given by 



H = {l + Xx^)?^ = ^[VT+X^P 



Contrarily to the naive expectation of ordering ambiguities, however, the 
usual procedure of canonical quantization does not present any ambiguity be- 
cause, as pointed out before, the linear operator (we put h=l) 

d 



P= + Ax2 ^ 

ox 

is self adjoint in the space X^(M, d/x), and therefore the quantum Hamiltonian 
operator turns out to be 

In presence of an interaction described by a potential Vi{x) things works 
similarly and the Hamiltonian is then 
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3 The factorization method 



Hereafter, as the configuration space is one-dimensional, we will use the notation 
d/ dx instead of the more traditional for the vector field d/dx. 

Let us try to determine a function W{x), called super-potential function, in 
such a way that the operator A and its adjoint operator A\ given by 

A =-L(vTW^ + H'(x)), 
are such that Hi = A^ A [13], i.e. 

Hl = A^A = l\-^/^ + Xx^ ■^ + w{x)]\^/^ + Xx^ ■^ + w{x) 

2 1 dx J L dx J 

Therefore, the super-potential function W must satisfy the following Riccati 
type differential equation 

Vl + W' -W"^ + 2Vi=0 . 

Once such a factorization is obtained, we can define a new quantum Hamil- 
tonian operator 



H2 = AA^ 



dx 



W{x) 



dx 



W{x) 



which is called the partner Hamiltonian. The new potential V2 is given in terms 
of the super-potential W by 



The operator A is such that A Hi = H2A while A^ is such that A^ H2 = 
HiA^. This shows that if |^) is an eigenvector of Hi corresponding to the 
eigenvalue E, then when A\^) 7^ 0, A\^) is an eigenvector of H2 corresponding 
to the same eigenvalue, because 

^2^1*) =^4^1 1*) =EA\^) , 

and similarly, when |$) is an eigenvector of H2 corresponding to the eigenvalue 
E and such that 74^]$) ^ 0, then |$) is an eigenvector of Hi corresponding 
to the same eigenvalue E. In other words, the spectra of Hi and H2 are almost 
identical, the only difference appearing when either \^) is an eigenvector of Hi 
but A \^) = 0, or 1$) is an eigenvector of H2 for which |$) = 0. 

As an important first remark, the potential function V is only defined up 
to addition of a constant, and therefore all the preceding expressions can be 
extended to the case in which H is replaced by + c, where c is a constant. As 
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another remark, it happens very often that some parameters may appear in the 
expression of the potential function, and therefore the super-potential function 
W will also depend of the values of such parameters. The most important case 
is when the explicit forms of the potential and its partner are quite similar and 
only differ in the values of the parameters, and then the problem is said to have 
shape invariance. As we will see this is the case for the quantum nonlinear 
oscillator we are considering in this paper. 



4 The quantum nonlinear Oscillator 



In the particular case of the nonlinear harmonic oscillator for which the Haniil- 
tonian is given by 



H 



1 



{1 + Xx')pi + 
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1 + Aa;2 



the quantum Hamiltonian operator will be 



dx"^ dx I + Xx^ ^ 

Now, if for any real number /3 we define the linear operator in (R, d^) 

f3x 



A = 



V2 V dx Vl + Ax2 J 



for which its adjoint operator is 



AT 



then, we find that 



(3x 



dx Vl + Ax2/ 



(4) 



(5) 



(6) 



1 + Aa;2 



2^ 



where by simple comparing with Q we conclude that the hamiltonian iJ( = 
Hi — (l/2)/3 admits a factorization 

H[ := At A (7) 
whenever the parameters a and /3 are related by = (3{I3 + A), i.e. 

/3 = (a- Va^TI^) . 



Note that 



lim B = a 



5 



Now for the partner hamiltonian iJj '■= AA^ we find 
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(8) 

Coming back to the general case, we recall that when a quantum Hamiltonian 
Hi[a) depending on a parameter a admits a factorization in such a way that 
the partner Hamiltonian H2 (a) is of the same form as Hi (a) but for a different 
value of the parameter a, it is usually said that there is shape invariance. More 
specifically, the quantum Hamiltonian Hi{a) admitting a factorization Hi{a) = 
A'' (a) A{a) is said to be shape invariant when there exists a function / such 
that for the partner H2(a) — A{a) A^{a) we have: 

H2{a) ^ Hiiai) + R{ai) , (9) 

where ai = f{a) and R{a) is a constant depending on the parameter a. In this 
case, it was shown in [14] and [15] that there is a method for exactly computing 
all the spectrum of Hi (see e.g. [9] for a modern approach). The bound state 
l^'o(a)) is found by solving A(a)|\['o(a)) — 0, and has a zero energy. Then, 
using © we can see that |^o(Q!i)) is an eigenstate of H2{a) with an energy 
El = R{ai), because 

i?2 (a) 1^-0(01)) - (Hiiai) + i?(ai))|*o(ai)) - i?K)|«'oK)) • (10) 

Then, A'l'(a)|^o(«i)) is the first excited state of Hi{a), with an energy 
El = i?(ai), because 

Hi{a)A^a)\^niai)) - A^ ia)H2{a)\^o{ai)) = A^(a)(i/i(ai) + i?(ai))|*o(ai)) 

- i?(ai)At(a)|^'o(ai)). 

This process should be iterated and we will find the sequence of energies for 
Hi{a) 

k 

£;fc=^i?(a,), £^0 = 0, (11) 
i=i 

the corresponding eigenfunctions being 

|*„(ao)) - At(ao)At(ai)...At(a„_i)|*o(«„)), (12) 

where = a and a^+i = f{<^j), namely, = f^{ao) — f^{ct)- 

Now we can apply this process to the case we were considering, the parameter 
being /3, because when comparing H'l given by 10 with its partner ((Sj) we see 
that as 



1 



2 



, d 1 1 _ . . . / \ 1 



H'i{f3-\)^-- + (/3_A)/? — — -3 --(/3-A), 



dx\ 2^' Kl + Xx"^ J 2 
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then 

and, therefore, 

where / is the function /(/?) = /? — A. If ii is the function defined by R{f3) = 
P + (1/2) A, then we see that 

H!,{P) = H[{0^) + R{f3i) . 

Therefore, as the quantum nonlinear oscillator we are considering has a shape 
invariance, we can develop the method sketched before for finding both the 
spectrum and the corresponding eigenvectors. 

The usual choice of parameters for these shape invariant systems is such that 
the function / corresponds to a displacement by one unit and then this suggests 
us to use the parameter 7 = /3/A instead of f3. 



5 The spectrum of the quantum nonlinear oscil- 
lator 

Our starting point should be the bomid state |'I'o(/3)) of the Hamiltonian H[. 
This eigenvector is determined by the condition A{P) |^o(/3)) = 0. More specif- 
ically, we should solve the differential equation 



ax 



1 + Aa;2 



*o = 



and therefore the wave function of the fundamental state must be proportional 
to 

The energies of the excited states will be 

E[ = R{0i) = /3 - A + A/2 = /3 - A/2 , 
and iterating the process we get 

K = f^i?(/3,) = i:(/3fc + ^)=E(/3-Afc + ^) , 

fe=l fe=l ^ ^ k=l ^ ^ 



and therefore, 



=n/? + A 



2 h J 



n/?-— A. 
2 



7 



The energy of the eigenstates of Hi = H{ + (l/2)/3 wiU be given by 



E^=np-—X+-p. 



This relation can also be written as 




) 



2 



/3(/?-A) 



+ 



2A 



The method developed in [14] and [15] also provides us the corresponding 
eigenfunctions as 



There is a clear difference between the cases A > and A < 0. In fact, let 
us first remark that the lowest value for E'^ is Eq = 0. Therefore, if A > only 
those values n are allowed for which 



On the contrary, when A < all natural numbers arc allowed for n. It should 
also be remarked that as long as X ^ 0, i.e, for both signs of A, the eigenvalues 
are not equally spaced. 

6 Final Comments and Outlook 

In this paper we have analyzed a quantum version of a nonlinear quantum 

oscillator with an amplitude dependent angular frequency. In the limit when 
the parameter A tends to zero we recover the usual harmonic oscillator. So, we 
first note that if A = 0, then a = jS, and if we take into account that 



wc sec that the ground state we have found coincides with that of the corre- 
sponding quantum linear harmonic oscillator. The limit, when A goes to 0, 
of operators A and and the eigenvalues of the Hamiltonian are the corre- 
sponding annihilation and creation operators and the energy eigenvalues for the 
harmonic oscillator. 

Another interesting issue concerns the behaviour of the position of the energy 
levels when the parameter A changes. Wc remarked that equispacing of the 
levels, a characteristic property of the usual harmonic oscillator, no longer holds 
when A 7^ 0. In particular, the fundamental level is given by Eq = 0/2; while 
a superficial reading would suggest this is a constant independent of A wc must 
stress this is not the (3 itself has been determined as a solution of the 

quadratic equation = /3(/J + A), whose appropiate determination is 



= ^t(/3)At(/3i) . ..A\p„_i)M0„) ■ 





= 



-A + VA^ + 4a2 
2 
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Thus, if we consider the nonlinear oscillator with fixed 'strcnght' a and allow the 
nonlinear parameter to vary, the fundamental level will actually depend on A. 
The fundamental level Eq (for n = 0) has the standard oscillator value a/2 for 
A = 0; and is a positive and always decreasing function of A. For large negative 
values of A it approaches the large positive asymptotic regime Eq{X) ~ —A/2, 
while for large positive values of A it tends to zero. 

It is a well-known fact that canonical transformation and quantization do 
not commute. We have made a canonical point transformation in such a way 
that the new coordinate is the one rectifying the vector field generating the 
isometrics of the metric associated to the Lagrangian describing the system and 
then the corresponding invariant measure is but the differential of such variable 
and the Hamiltonian becomes the square of the momentum, eliminating then 
order ambiguities. With this quantization procedure wc have found that the 
spectrum and the corresponding eigenvectors can be easily found using the fact 
that the Hamiltonian admits such factorization that it is shape-invariant. 
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